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ABSTRACT 

5 This paper presents a simple analytical framework for the dynamic response of cirrus to a 

6 local radiative flux convergence, expressible in terms of three independent modes of cloud 

7 evolution. Horizontally narrow and tenuous clouds within a stable environment adjust to 

8 radiative heating by ascending gradually across isentropes while spreading sufficiently fast 

9 SO as to keep isentropic surfaces nearly flat. More optically dense clouds experience very 

10 concentrated heating, and if they are also very broad, they develop a convecting mixed layer. 

11 Along isentropic spreading still occurs, but in the form of turbulent density currents rather 

12 than laminar flows. A third adjustment mode relates to evaporation, which erodes cloudy 

13 air as it lofts. The dominant mode is determined from two dimensionless numbers, whose 

14 predictive power is shown in comparisons with high resolution numerical cloud simulations. 

15 The power and simplicity of the approach hints that fast, sub-grid scale radiative-dynamic 

16 atmospheric interactions might be efficiently parameterized within slower, coarse-grid climate 

17 models. 
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18 1. Introduction 



19 Cloud-climate feedbacks remain a primary source of uncertainty in climate forecasts 



(Dufresne and Bony 2008), mainly because clouds both drive and respond to the general 
circulation, the hydrological cycle, and the atmospheric radiation budget. Unlike fields of 
water vapor, clouds evolve quickly, so their radiative forcing and dynamic evolution are 
highly coupled on time and spatial scales that cannot be easily resolved within Global Cli- 
mate Models (GCMs). For faithful reproduction of large-scale climate features, resolving 
radiatively driven motions on sub-grid scales may be at least as important as accurately 



25 representing mean grid-scale fluxes (Cole et al. 2005). 



Radiative flux convergence and divergence within cloudy air is normally thought to pro- 



28 duce vertical lifting and mixing motions (Danielsen 1982; Ackerman et al. 1988 Lilly 1988 



Jensen et al. 1996 Dobbie and Jonas 2001). What is often overlooked is that clouds with a 



finite width also adjust to radiative heating by spreading horizontally, especially if the heat- 



si ing is concentrated in a thin layer at the cloud top or bottom (Garrett et al. 2005, 2006). 

32 Such radiatively driven mesoscale circulations have been identified within thin tropopause 

33 cirrus, and they are thought to play a role in determining the heating rate of the upper 



troposphere (Durran et al. 2009) and in stratospheric dehydration mechanisms (Dinh et al. 
2010). Jensen et al. (2011) suggest that radiative cooling can help to initiate thin tropopause 
cirrus formation, while subsequent radiative heating in an environment of weak stability can 
induce the small-scale convection currents that are required to maintain the cloud against 
gravitational sedimentation and vertical wind shear. 

Where these recent studies directly simulated the highly interactive and complex nature 
of cloud processes, an alternative and perhaps more general approach is to start with simple, 
analytical and highly idealized models that emphasize specific aspects of the relevant physics. 
Here, we look at the respose of cirrus clouds to local thermal radiative flux divergence within 
cloud condensate. The discussion that follows largely neglects precipitation, synoptic scale 
motions, and shear dynamics to facilitate description of a simple theoretical framework within 



45 a parameter space of two dimensionless numbers. A similar approach has been employed 

46 previously to constrain small-scale interactions between diabatic heating and atmospheric 



dynamics (Raymond and Rotunno 1989), including situations where radiation is absorbed 



48 by horizontally infinite clouds (Dobbie and Jonas 2001). Here, we extend consideration to 

49 radiatively absorptive layers that have finite lateral dimensions, an ingredient that turns out 

50 to be critical for predicting the evolution of cloud size and cross-isentropic motions. The 

51 broad intent of this study is to provide insight into how clouds respond to rapid, small-scale 

52 radiative heating in a way that might be better parameterized within large scale, coarse-grid 

53 models such as GCMs. 



54 2. Non-equilibrium radiative- dynamic interactions in cir- 

55 rus 

56 The starting point is to consider a microphysically uniform, optically opaque cloud that 

57 is initially at rest with respect to its surrounding, characterized by a stably stratified at- 

58 mosphere with a virtual potential temperature 9^ that increases monotonically with height 

59 (Figure [T]). The arguments described below apply equally to cloud base and cloud top, dif- 

60 fering only in sign of forcing. However, for the sake of simplicity the focus here is on cloud 

61 base. 

62 At cloud base, cloudy air has a lower brightness temperature than the brightness temper- 

63 ature of the ground and lower tropospheric air that is below it. This radiative temperature 

64 difference drives a net flow of radiative energy into the colder cloud base, effectively due to 

65 a gradient in photon pressure, that can be approximated as 

AF^ei - 4af ^Af (1) 

66 where a is the Stephan-Boltzmann constant, Tc is the cloud temperature, and AT is the 

67 effective brightness temperature difference between the lower tropospheric air and cloud 
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68 base. 

69 Provided the cloud is sufficiently opaque to act as a blackbody, radiative energy is de- 

70 posited within a layer of characteristic depth h at the base of the cloud that is smaller than 

71 the depth of the cloud itself. The magnitude of h can be obtained by considering that the 

72 thermal emissivity is given by 

£ ~ 1 - expi-Tabs/li) (2) 

73 where Tabs is the absorption optical depth and /i is the quadrature cosine for estimating the 



74 integrated contribution of isotropic radiation to vertical fluxes. Usually /i ~ 0.6 (Herman 



75 1980). The absorption optical depth is determined by the cloud ice mixing ratio qi, as well 

76 as the ice crystal effective radius Vg, through Tabs = k{re)qipAz where k is the mass specific 

77 absorption cross-section density, p is the density of air, and Az is the vertical path length 

78 through which the radiation is absorbed. The depth h is the e-folding path length for the 

79 attenuation such that Tabs/p = 1: 

k{re)qiP 

so Assuming an effective radius of 20 pm, the value for k{rg) in cirrus is approximately 0.045 



81 m^g ^ ( KnoUenberg et al. 1993). Taking, for example, qi values of 1 g kg ^ that have been 



82 observed in medium sized cirrus anvils in Florida ( Garrett et al.||2005 ), the depth h would be 



83 about 30 m. As a contrasting example, a cloud with q^ values of 0.01 g kg ^, similar to those 



84 observed in thin cirrus (Haladay and Stephens 2009), would have a radiative penetration 



85 depth h of about 3000 m. Thus, the deposition of radiative enthalpy in this layer increases 

86 its temperature at rate 

n = = — ~ ^ = AaT^^^^AT (4) 

dt pcp dz pcph pcp 

87 where Cp is the specific heat of the air. To ffist order, heating rates are proportional to the 

88 radiative temperature contrast and the cloud ice mixing ratio. 
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89 a. Dynamic Adjustment to Diabatic Heating 

90 The total flow of upwelling radiative energy into a cloud is proportional to AFnet and the 

91 normal cloud horizontal cross-section, which is of order where L is the cloud horizontal 

92 width. Deflning the initial, neutrally buoyant, ground-state for the gravitational potential 

93 energy density of the cloudy air within the volume hL'^ as = E^q/ [hL"^) (Figure [T]), then 

94 an accumulated flow of energy into the volume increases the gravitational potential energy 

95 density to + Ae at rate dAe/dt = R AFnet / (cph) , where R is the gas constant for air. The 

96 remaining fraction (cp — R)/cp = c^/cp of the radiative enthalpy deposited in the cloud goes 

97 towards increasing the rotational and translational energy density of the cloudy air within 

98 the layer. Conceptually, it is useful to consider the gravitational increase as an increase in 

99 the pressure gradient that is available to drive fluid dynamic motions: pressure gradients 

100 have units of energy density. 

101 The increase in the potential energy density within the volume hL"^ allows work to be 

102 done against the overlying gravitational static stability to create a mixed-layer with, on 

103 average, near constant 6i,. Thus any newly absorbed thermal energy becomes redistributed 

104 through a mixed layer depth 6z that is larger than the radiatively absorbing layer of depth 

105 h (Figure [T]). This is important, because it has the effect of diluting the density of newly 

106 added radiative energy through a factor of 6z/h such that: 

dAe ^ RAFnet 
dt Cp6z 

107 As required by the second law of thermodynamics, equilibrium is restored through relaxation 

108 of the buoyant potential energy density perturbation Ae to zero, leading to kinematic flows 

109 (Figure [2]). 

no There are two basic modes for relaxation of the buoyant energy density. The available 

111 gravitational potential energy density Ae can be expressed as the density p of the air at a 

112 given buoyant potential energy density, multiplied by the buoyant potential per unit mass 
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113 of air that is available to drive flows 

~ N'^5z^ (6) 

114 Here, N is the buoyancy frequency, which is related to the local stratification through 

e,dz ^ ' 

115 It follows that 

Ae = pN'^5z^ (8) 

116 Dynamic relaxation of the radiatively induced perturbation Ae can proceed in either of 

117 two ways. At constant density, the heated volume L'^5z can be raised to higher gravitational 

118 potential. Alternatively, the air expands outwards along a constant potential surface. 

(ilnAe 91nAE d\np 

119 Given Eq. |6] and that 

mm 



where m/5z is fixed (i.e., no entrainment of mass across the mixed-layer boundary), Eq. [9| 
can be rewritten as 

(ilnAe d\n5z dhiL 

= '^^ ^^^^ 

=a5z - OIL (12) 
120 where a^z and represent instantaneous rates of adjustment. 



121 Eq. Y\_ has several implications. The buoyant potential energy density Ae within the 

122 mixed-layer volume L'^5z can grow due to the continuing radiative fiux deposition within the 



123 volume (the positive first term in Eq. 11). Or, it can decay through horizontal expansion 



(the negative second term in Eq. 11). In the first case, if the width of the cloud L is held 



125 constant, the mixed-layer deepens into stratified cloudy air above it at rate 

d{5z) . dejdt Ugh 



dt dejdz e^N'^dz 
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(13) 



126 where the factor of h/5z arises from the dilution of potential energy through a depth larger 

127 than the absorptive layer where initially, 5z = h. Mixed-layer growth rates slow with time. 



128 The solution to Eq. [13] as a function of time At is 

^^-i^M" (14) 

129 Alternatively, if the depth of the mixed-layer Sz is fixed, then the potential energy density 

130 relaxes towards equilibrium by smoothing out horizontal pressure gradients between the 

131 cloudy mixed-layer and clear sky beside it. It does this through expansion of the volume 

132 L'^6z along constant potential surfaces (or isentropes) into the lower potential energy density 

133 environment that surrounds the cloud. This density current outflow occurs at speed 

dL 

Umix = -g^ \sz^ N6z (15) 

134 which results from the conversion of the gravitational potential energy of order N'^Sz'^ into 

135 kinetic energy of order M^j^. 

136 To assess the relative importance of cross-isentropic adjustment to along-isentropic spread- 

137 ing, a dimensionless number 5* can be defined as the ratio of the two rates asz and in Eq. 

138 [TT] From Eq. [13} radiative heating increases the mixed-layer gravitational potential energy 

139 density at rate 

dlnSz 2'Hgh 



140 From Eq. 15, the rate of loss of potential energy density due to expansion of the mixed-layer 

141 laterally into the clear-sky surroundings is 

91nL , N5z 
a, ~ 2— |,= 2— (17) 

142 For a cloud that is initially at rest, in which case a mixed layer has not yet developed, 

143 then radiation deposition remains concentrated within the layer 6z h and from Eqs. [16] 



144 and 17, the ratio of these two rates can be defined by a dimensionless "Spreading Number'' 

asz d\ia6z/dt\L HgL 



S 
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145 If S* > 1, then the potential energy density within the layer L'^5z increases due to radiative 

146 flux deposition at a rate asz that is faster than the rate at which gravitational relaxation 

147 can reduce the disequilibrium in potential energy density through horizontal flows into sur- 

148 rounding clear air. Isentropic surfaces at cloud base cannot stay flat, but rather are deformed 

149 downward by the radiative heating. This deformation creates a deepening turbulent mixed 

150 layer that gradually grows into the overlying static stability of the atmosphere as the square 



161 root of time (Eq. 14). Meanwhile, the mixed- layer spreads outward along isentropic surfaces 



152 at rate u^ix (Eq- 15). 



153 By contrast, when < 1, adjustment through isentropic spreading is sufficiently rapid 

154 that isentropic surfaces stay approximately flat. Cloud motions stay laminar rather than 

155 becoming turbulent. The mixed-layer horizontal expansion given by (Eq. [I7| ) decreases 

156 the potential energy density faster than the rate asz (Eq. [16]) at which potential energy 

157 density is deposited at cloud base through radiative flux convergence. The potential energy 

158 density at cloud base does not increase and does not overcome the overlying static stability. 

159 Rather the cloud simply lofts across isentropic surfaces at speed 

_ n _ ng 

160 Dimensional continuity arguments require that the cloud spreads laterally along isentropes 

161 at speed 

L UgL . , 

162 h. Evaporative Adjustment 

163 The above describes two modes for how radiative flux deposition can create pressure, 

164 or potential energy density gradients that drive cloud-scale motions. A third possibility 

165 for adjustment is that local radiative heating may result in microphysical changes where 

166 temperature is maintained, but condensate evaporates or condenses. 

167 Assuming that all absorbed radiative energy goes towards evaporation at cloud base. 



168 and that there is no lag associated with the diffusion of vapor away from ice crystals, then 

169 ice evaporates at rate pLgdqi/dt = AFnet/h, where Lg is the latent heat of sublimation. 

170 Substituting Eq. [3] for h, radiative heating evaporates cloud base at rate 



dlnqi 
dt 



(21) 



171 Note that if there were net radiative flux divergence, as might be expected at the top of a 

172 thermally opaque cloud, then net cooling would lead to condensation. 



173 The ratio of Oevap (Eq- 21 ) to asz (Eq- 16 ) implies a dimensionless "Evaporation Number 



174 comparing the evaporation rate to the rate of laminar adjustment through cross isentropic 

175 ascent. In the initial stages of development, where 5z = h, 



176 or, substituting Eq. [7] and Eq. |4] for the heating rate T-L 

ijsQi az 

177 The susceptibility to evaporation depends only on the cloud microphysics and the local 

178 static stability, and not, in fact, on the magnitude of the heating. Provided E > 1, cloud 

179 base evaporates rather than lofts. However, for values of -E < 1, cloud ascends faster than it 

180 evaporates and condensate is maintained. It is important to note here that the "Evaporation 

181 Number" E should only be considered if the "Spreading Number" S has values smaller than 

182 unity. If 5 > 1, the relevance of evaporation is less clear because a convective mixed-layer 

183 develops, in which case one would expect instead continual reformation and evaporation of 

184 cloud condensate as part of localized circulations within the mixed-layer. The more relevant 

185 comparison might be to rates of turbulent entrainment and mixing. 
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186 3. Numerical Model 



187 To test the suitabihty of the dimensionless "Spreading" and "Evaporation" numbers S 

188 and E for determining the cloud evolutionary response to local diabatic heating, we made 

189 comparisons to cloud simulations from the University of Utah Large Eddy Simulation Model 

190 (UU LESM) ( Zulauf||2001 )■ An LES model is used because the resolved scales are sufficiently 



191 small to represent turbulent motions, convection, entrainment and mixing, and laminar flows. 

192 

193 The UU LESM is based on a set of fully prognostic 3D non-hydrostatic primitive equa- 



194 tions that use the quasi-compressible approximation (Zulauf 2001). The model domain was 

195 placed at the equator, = 0°, to eliminate any Coriolis effects. Even in the largest domain 

196 simulations, the maximum departure from the equator (50 km) is sufficiently small as to 

197 justify not including the Coriolis effect in the model calculations. 

198 The horizontal extent of the domain was chosen to contain the initialized cloud as well 

199 as to allow sufficient space for spreading of the cloud during the model run. The UU LESM 

200 model employs periodic boundary conditions such that fluxes through one side of the domain 

201 (moisture, cloud ice, turbulent fluxes, etc.) enter back into the model domain from the 

202 opposite side. Here, the horizontal domain size is case dependent but chosen to be sufficiently 

203 large as to minimize "wrap around" effects. Horizontal grid size was chosen to be 30 m to 

204 match the minimum value for vertical penetration depth of radiation into the cloud, but it 

205 increased to 100 m for cases that required particularly large and computationally expensive 

206 domains. 

207 The vertical domain spanned 17 km and included a stretched grid spacing. The highest 

208 resolution for the stretched grid was placed at the center of the initial cloud with grid size 

209 of 30 m. The vertical resolution decreased logarithmically to a maximum grid spacing of 

210 approximately 300 m at the top of the model and approximately 400 m at the surface. A 

211 sponge layer was placed above 14 km to dampen vertical motions at the top of the model 

212 and to prevent reflection of gravity waves off the top of the model domain. The model time 
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213 step for dynamics was between 1.0 and 10.0 s and was chosen to be the largest time step 

214 that was computationally stable. 

215 For radiative transfer, the UU LESM model uses a plane parallel broadband approach, 



216 using a 5- four stream scheme for parameterization of radiative transfer (Liou et al. 1988), 



217 based on the correlated k-distribution method (Fu and Liou 1992). Radiative transfer cal- 

218 culations were performed at a time step of 60 s. Only thermal radiation was considered in 

219 this study. 

220 For all cases examined, the model was initialized with a standard tropical profile of 

221 temperature and atmospheric gases with a buoyancy frequency of approximately 0.01 s"^ 

222 throughout the depth of the model domain. Relative humidity was set in two independent 

223 layers. In the bottom layer of the model, which extends from the surface to 7.8 km, the 

224 relative humidity was set to a constant 70% with respect to liquid water. In the upper layer 

225 of the model, from 7.8 km upwards, which contained the cloud between 8.8 km and 11.3 km, 

226 relative humidity with respect to ice was set to a constant value of 70%. All clouds were 

227 initialized as homogeneous cylindrical ice clouds, as shown in Figure |3| Ice particles within 

228 the cloud were of uniform size with a fixed effective radius of 20 /xm and an initially uniform 

229 mixing ratio as prescribed by the particular case. Cloud radius was prescribed according 

230 to the particular case, but in each case the thickness was set to 2500 m with the cloud 

231 base set at 8.8 km. Cloud base was chosen such that the cloud top would be placed at 

232 approximately 200 mb, in rough accordance with the average cirrus anvil height indicated 



233 by the Fixed Anvil Temperature hypothesis (Hartmann and Larson 2002). Both the cloud 

234 and surrounding atmosphere were initialized to be at rest. No precipitation was allowed in 

235 any of the model simulations. Cloud particle fall speed was also neglected. All cases were 

236 run for one hour of model simulation time. 

237 Two cloud parameters were varied through several orders of magnitude in order to explore 

238 a wide parameter space of possible evolutionary behaviors. Cloud radius L was chosen to be 

239 100 m, 1 km, or 10 km. The ice-water mixing ratio qi was set to 0.01, 0.1, or 1 g kg~^. This 
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240 provided 9 unique combinations of cloud size and density, as described in Tables [T] and [2} 

241 that spanned a range of values of 5* and E, and included combinations that are sufficiently 

242 unstable that they are not observed naturally. 

243 Figure |4] shows the initial heating rate profiles for each value of qi used in this study 



244 calculated using the Fu and Liou ( 1992 ) radiative transfer parameterization. Note that the 

245 heating is confined to a narrower layer at cloud base as the ice water mixing ratio increases 

246 (Eq. The heating profiles for both the Qi = 0.01 g kg~^ and = 0.1 g kg~^ cases closely 



247 match the calculated heating rate profiles from Lilly (1988 ). However, the heating rate profile 

248 for the = 1 g kg^^ case, which Lilly did not model, shows an order of magnitude increase 

249 in the heating and cooling rates to several hundred K day~^, confined almost exclusively to 

250 the top and bottom of the cloud, with virtually no heating in the interior. 

251 For cases with qi = 0.01 g kg~^, the radiative penetration depth h is 3300 m, which is 

252 deeper then the 2500 m cloud depth. However, the heating profile is nearly linear through 

253 the depth of the cloud with heating at cloud base and cooling at cloud top. Thus, in cases 

254 where qi is 0.01 g kg~^, the radiative penetration depth h is assumed to be half the cloud 

255 depth, or 1250 m, for the purposes of calculating S and E. 



256 4. Results 

257 In the parameter space of S and E described by Tables [T] and |2| tenuous and narrow 

258 clouds with low values of ice water mixing ratio qi and cloud width L have values of the 

259 spreading number S that are less than 1. Theoretically, such clouds are expected to undergo 

260 laminar lifting and spreading. Tenuous clouds with large values of E and small value of 5* 

261 are expected to evaporate at cloud base. Optically dense and broad clouds with large values 

262 of qi and L have values of S much larger than 1, and are expected to favor the concentration 

263 of potential energy density in a thin layer at cloud base, leading to turbulent mixing and 

264 erosion of stratified air within the cloudy interior. 
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265 In what follows, numerical simulations are performed to test the validity of the dimen- 

266 sionless numbers 5* and E for predicting cloud evolution. Cases that describe the parameter 

267 space in S will be discussed first, since values of E are relevant only for scenarios with 

268 5 < 1 where mixed-layer development is not the primary response to local diabatic radiative 

269 heating. 



270 a. Isentropic Adjustment 

271 Simulations of clouds with values of S* < 1 are expected to show cross-isentropic ascent 

272 of cloud base in response to local diabatic radiative heating and, through continuity, laminar 

273 spreading. Effectively, the loss of potential energy out the sides of the cloud (due to material 

274 flows) is sufficiently rapid to maintain nearly flat isentropic surfaces within the original 

275 cloud volume. Equivalently, cross-isentropic ascent is sufficiently slow that the consequent 

276 horizontal pressure gradients can be equilibrated through laminar spreading while keeping 



277 isentropic surfaces approximately flat (Eq. 17). 



278 A good example of this behavior is shown in a simulation of a cloud with L = 1 km and 

279 Qi = 0.1 g kg~^. This case has a value oi S = 0.033, which implies that the primary response 

280 to radiative heating should be adjustment through ascent across isentropic surfaces. Figure 

281 |5] shows the isentropes, or contours in 6^. The isentropes remain approximately flat and 

282 unchanged from their initial state in response to the cross-isentropic flow of cloudy air. As 

283 shown in Figure [6} the simulated cloud undergoes rising at cloud base and sinking at cloud 

284 top, while spreading horizontally. 



b. Mixing 

Clouds with values of S* > 1 are not expected to be associated with laminar motions. 
Instead, radiative heating bends down isentropic surfaces so rapidly as to create a local 



instability that cannot be restored sufficiently rapidly by laminar cloud outflows (Eq. 16). 
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289 Radiative heating is sufficiently concentrated to initiate turbulent mixing that produces a 

290 growing mixed-layer. Unlike the 5* < 1 case, isentropes do not stay flat. 

291 An example, shown in Figure [7| is for a simulated cloud that has initial condition values 

292 of cloud radius L = 10 km and ice water mixing ratio = 1 g kg~^. Since S = 1300, it 

293 is expected that the potential energy density at cloud base will increase at a rate that is 



294 faster than the loss rate of potential energy through cloud lateral expansion (Eqs. 16 and 



17). A mixed layer will develop because the deposition of radiative energy creates buoyancy 



296 that does work to overcome the static stability of overlying cloudy air and create a mixed 



297 layer. Meanwhile the mixed-layer expands with speed Umix = N6z (Eq. 15), where 6z is the 

298 mixed-layer depth and is the static stability of air surrounding the cloud. 

299 The numerical simulations reproduce these features. A mixed-layer can be seen in the 

300 9y profile plotted in Figure [7} showing the average cloud properties after 1 hour of model 

301 simulation. This proffie is a horizontally averaged profile taken within 9 km of cloud center. 

302 On average, the mixed-layer exhibits a nearly adiabatic profile in At 1 h simulation time, 

303 the mixed-layer at cloud base is nearly 800 m deep. The mixed layer expands horizontally 

304 along isentropes, as seen in Figure |8| The "bowl" shaped spreading of the cloud is because 

305 intense radiative heating at cloud base bends isentropic surfaces downward. 

306 This mixed-layer development and spreading can also be seen in cross sectional plots of 

307 qi in Figure |9] There is a mixed-layer at both cloud base and cloud top with darker shading 

308 indicating where drier air has been entrained from below or above. Note that cloud base 

309 and cloud top remain at roughly constant elevation. In Figure [5} for a case where 5* ^ 1, 

310 radiative flux convergence at cloud base drives cross-isentropic laminar ascent. In this case, 

311 where 5* ^ 1, laminar ascent does not occur. Instead, cloud base remains nearly at its 

312 initial vertical level and there is formation of a turbulent mixed layer that spreads outward 

313 along isentropes. Notably, the mixed-layer circulations at cloud base have a mammatus-like 



314 quality to them, something we have discussed more extensively in Garrett et al. (2010). 

315 These behaviors can be quantified by examination of the rapidity of development of a well- 
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316 mixed layer at cloud base. If the dominant mode of evolution is cross-isentropic lofting, then 

317 vertical potential temperature gradients should remain relatively undisturbed. Conversely, 

318 if mixing is the dominant response, then potential temperature will evolve to become more 

319 constant with height. 

320 Table |3] shows the cloud domain-averaged, logarithmic rate of decrease in the static 

321 stability din N"^ /dt, where A^^ oc dO^/dz. Calculations are evaluated for the lowermost 80 

322 m of the cloud within the initial 360 s of simulation time. The destabilization of cloud 

323 base reflects the magnitude of the Spreading Number S (Table [T]), with large values of S 

324 demonstrating the most rapid rates of mixed-layer development. 



325 c. Evaporation 

326 Cloud bases with S* < 1 and E > 1 are expected to evaporate more quickly than they loft 

327 across isentropes (Table [2|. For example, for a cloud with L = 1 km and qi = 0.01 g kg"^ , 

328 the calculated value of the Spreading Number S is 0.0011, and the value of the Evaporation 

329 number E is 150. Based on these values, the expected evolution of cloud base would be 

330 gradual evaporative erosion of cloud base. 

331 To quantify the importance of evaporation to cloud evolution, the rate of change in cloud 

332 mass d\nm/dt, where m is the mass of cloud ice, was calculated over the first 180 s of 

333 simulation, but only within the lower layer in which radiation from the surface is absorbed, 

334 h, rather than the entire cloud. The absorptive layers were taken to be 30 m, 300 m, and 

335 1250 m for cloud ice water mixing ratios of 1 g kg~^, 0.1 g kg~^, and 0.01 g kg~^ respectively. 



336 From Eq. 21 for aevap, the anticipated evaporation rate at cloud base is approximately 

337 7 h^^ based on the modeled net flux absorption AE of 74 W m~^ within the absorption 

338 layer h. Table |4] shows maximum modeled evaporation rates that are nearly as large, at 

339 least where E is maximized and the cloud is narrow. However, rates of evaporation decrease 

340 with increasing cloud width L, perhaps because S increases and stronger dynamic motions at 

341 cloud base replace evaporated cloud condensate with newly formed cloud matter. In general, 
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342 however, tenuous cirrus clouds are most susceptible to erosion by evaporation at cloud base, 

343 particularly if they are not very broad. 



344 d. Precipitation 

345 While the role of precipitation has been excluded from these simulations in order to clarify 

346 the physical behavior, certainly natural clouds can have significant precipitation rates. An 

347 estimate of the relative importance of precipitation is briefly discussed here. 

348 The characteristic precipitation timescale apredp depends on the rate of depletion of 

349 cloud water by precipitation P and the average ice water content IWC. For example, in 

350 a cirrus anvil in Florida measured by aircraft during the CRYSTAL-FACE field campaign, 

351 the measured value of P was 0.05 g m~^ h~^ compared to values of IWC of 0.3 g m~^ 



352 (Garrett et al. 2005), implying a precipitation depletion rate Opredp = P/IWC = 0.15 h . 



353 For comparison, corresponding values for the radiative adjustment rates are ai — O-ll ^ 



354 (Eq. 17) and asz — 144 h~^ (Eq. 16). While development of a turbulent mixed layer is 

355 the fastest process, precipitation depletes cloud condensate at a rate that is comparable to 

356 ai, the rate at which gravitational equilibrium is restored through cross-isentropic flows and 

357 laminar spreading. 



358 5. Discussion 

359 We have separated the evolutionary response of clouds to local diabatic heating into 

360 distinct modes of cross-isentropic lifting, along-isentropic spreading, and evaporation of cloud 

361 condensate. A straightforward method has been described for determining how a cloud will 

362 evolve based on ratios of the associated rates. The dominant modes of evolution are outlined 

363 in Table |5l 

364 For example, cirrus anvils begin their life cycle as dense cloud from convective towers 

365 that have reached their level of neutral buoyancy ( Scorer|[l963 Jones et al.||1986 Toon et al. 
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2010). Such broad optically thick clouds are associated with high values of the spreading 



number S due to their large horizontal extent and high concentrations of cloud ice. Radiative 
flux convergence is confined to a thin layer at cloud base. Heating is so intense, and the cloud 
is so broad, that the cloudy heated air cannot easily escape by spreading into surrounding 
clear air. Instead, large values of S favor the development of a deepening mixed-layer. The 
mixed-layer still spreads, but in the form of turbulent density currents rather than laminar 
motions. 

However, as the cloud spreads and thins, the value of the spreading number S evolves. 
S is proportional to the heating rate "H, cloud width L, and inversely proportional to the 



square of the depth of the mixed layer 62"^ (Eq. 18). Cloud spreading increases the value 
of L, and this acts as a positive feedback on S. But as the cloud spreads, the mixed-layer 



depth increases as 

(Eq. h4, progressively diluting the impact of radiative heating on 
dynamic development by a factor of 6z/h. Thus, while cloud spreading increases S, this is 



379 offset by increasingly diluted heating rates within the mixed-layer (Eq. 18). 



From Eq. 18, S can be rewritten as 



S = A 



6z^ 



(24) 



where A = Tigh/O^N^ is assumed to be constant, assuming here that qi is fixed. Thus, the 
rate of change in 5* is given by 

dlnS , 



dt 



dlnL dlnSz 
-3- 



dt 



dt 



(25) 



383 From Eq. 13, and since dL/dt = Umix ~ N6z (Eq. 15), Eq. 25 can be rewritten as 

d\nS , N6z SNA 



dt 



6z^ 



384 Finally, from Eq. 14, if the mixed layer depth evolves over time as 5z = (NAtY^"^, Eq. 

385 becomes 



(26) 
26 



d\nS _ {N^Aty/^ 3 

dt Ln t 



(27) 
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386 Thus, the evolution of S is controlled by two terms, the first being a positive feedback 

387 related to cloud spreading, and the second being a negative feedback related to mixed-layer 

388 deepening. Provided that 

/ 9^2 \ 1/3 

t < t^a. = j (28) 

389 the negative feedback dominates, so that to a good approximation 

dlnS , . 

— — ~ -3 29 
dint ^ ^ 

390 which can be solved for the general solution 

(30) 

391 For a thick cirrus anvil with initial values of of 1 g kg~^, L of 10 km, and S of 1300, the 

392 value of A is 3510 m^ and tmax — 10 h. By comparison, from Eq. 30, the value of S rapidly 



393 drops to a value of approximately unity within time t ~ lOto- While the value of to is not 

394 explicitly defined, assuming that it is one buoyancy period 2n/N, then the time scale for 

395 the cirrus anvil to shift from turbulent mixing to isentropic adjustment is of order one hour. 

396 Because this timescale is much less than tmax, the anvil never manages to enter a regime of 

397 runaway mixed-layer deepening where Eq. |27] is positive. What is interesting is that this 

398 timescale for a convecting anvil to move into a laminar fiow regime is comparable to the few 



399 hours lifetime of tropical cirrus associated with deep-convective cloud systems (Mace et al. 



2006). A transition to laminar behavior seems inevitable. 



401 Figure 10 shows numerical simulations for the time evolution of S within the cloud 

402 base domain. These reproduce the theoretically anticipated decay at a rate close to the 

403 anticipated power law. The decay in 5* is dominated by mixed-layer deepening, which 

404 roughly follows the anticipated t^/^ power law. It may seem counter-intuitive, but it is 

405 deepening of a turbulent mixed-layer that allows for a transition to laminar behavior: current 

406 radiative fiux deposition becomes increasingly diluted in past deposition. Once an anvil 

407 reaches S ^ 1, the rate at which the mixed layer deepens becomes roughly equal to the 
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rate at which laminar flow restores gravitational equilibrium through spreading. At this 
point, the dynamic evolution of the cirrus anvil enters a new regime where it adjusts to 
any radiatively induced gravitational disequilibrium through either cross-isentropic lofting 



411 (Danielsen 1982 Ackerman et al. 1988) or evaporation (Jensen et al. 1996) 



As a contrasting example, contrail formations are typically optically thin and horizontally 
narrow. In some cases they can evolve into broad swaths of cirrus that persist for up to 17 
hours after initial formation and radiatively warm the surface ( Burkhardt and Krcher||2011 ). 
While we did not specifically model contrails in this study, the theoretical principles that we 
discussed can provide guidance for how they might be expected to evolve. 

Immediately following ejection from a jet engine, the contrail air has water contents of 



a few tenths of a gram per meter cubed (Spinhirne et al. 1998), contained within a very 



419 narrow horizontal domain (Voigt et al. 2010). In this case, the cloud can be characterized 



in an idealized sense by qi =1 g kg^^ and L =100 m (Table [T|. Since the expressions for 
A and tmax discussed above do not depend on the horizontal extent of the cloud L, their 
values are identical to those of the idealized anvil that was explored. However, the initial 
value of S does depend on L, and with an initial value of 13 it is one hundred times smaller 
than for the anvil case. Since the initial value for S is still larger than unity, it should be 
expected that the contrail cirrus will be able to sustain radiatively driven turbulent mixing 



426 in its initial stages. However, from Eq. 30 , 5* should be expected to decline to unity in about 



20 minutes, at which point more laminar circulations take over that allow for the contrail 
cloud to spread laterally while slowly lofting across isentropes. 



429 6. Conclusions 

430 In this study, the evolutionary behavior of idealized clouds in response to local diabatic 

431 heating was estimated from simple theoretical arguments and then compared to high res- 

432 olution numerical simulations. Simulated clouds were found to evolve in a manner that 
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433 was consistent with expected behaviors. Dense, broad clouds had high initial values of a 



spreading number 5* (Eq. 18) and formed deepening convective mixed-layers at cloud base 



that spread in turbulent density currents. The mixed-layers were created because isentropic 
surfaces were bent downward by radiative flux convergence to create a layer of instability. 



The mixed-layer deepened at a rate a^^ (Eq- 16 ) that was much faster than the rate at which 
the potential instability could be restored through along-isentropic outflow into surrounding 



439 clear (Eq. For particularly high values of S, the mixed-layer production 



ifrom radiative heating was so strong as to create mammatus clouds at cloud base (Garrett 



441 et al. 2010). 



Tenuous and narrow clouds with initial values of < 1 displayed gradual laminar ascent 
of cloud base across isentropic surfaces while the cloud spread through continuity into sur- 
rounding clear sky. Isentropic surfaces stayed roughly fiat because the rate of along-isentropic 



spreading a^, (Eq. 17) was sufficiently rapid compared to the rate of cross-isentropic lifting 



af,z (Eq. 16 ) that isentropic surfaces in the cloud were continuously returned to their origi- 
nal equilibrium heights. Clouds with low values of S and also high values of an evaporation 



448 number E (Eq. 23 ) tended to evaporate quickly because the rate at which cloud condensate 



evaporated agt^ap (Eq- 21) was much faster than the rate at which the cloud layer rose in 



cross-isentropic laminar ascent ai^ (Eq. 17) 



For clouds with values of S that are initially high, the tendency is that S falls with time as 
the convergence of radiative flows at cloud base becomes increasingly diluted in a deepening 
mixed-layer. We found that dense cirrus anvils with a large horizontal extent remain in a 
mixed-layer deepening regime for nearly an hour before shifting across the 5 = 1 threshold 
into a cross-isentropic laminar lofting regime. Contrail cirrus are expected to make the same 
transition, but in a matter of tens of minutes. 

It is important to note that the precision of any of these results is limited by the sim- 
pliflcations that were taken. Most important is that no precipitation was included in the 
numerical simulations, so simulated clouds presumably persisted longer than if precipitation 
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460 were included. Also, single-sized ice particles were used rather than a distribution of ice par- 

461 tide sizes. Gravitational sorting would result in a higher concentration of larger ice particles 



462 ipear cloud base and a higher concentration of small ice particles near cloud top (Garrett 



et al. 


2005 


Jensen et al. 


2010 



464 Nonetheless, it has been shown that local diabatic heating heating can drive dynamic 

465 motions and microphysical changes that are at least as important as precipitation, and easily 

466 predicted from the simple calculation of two dimensionless numbers. A practical future 

467 application of this work might be improved constraints of the fast, smale-scale evolution of 

468 fractional cloud coverage within a GCM gridbox, limiting the need for explicit, and expensive, 

469 fluid simulations of sub-grid scale processes. 
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543 List of Tables 

1 Spreading number S=^ = EH 

545 2 Evaporation number E = = EH 

546 3 Rate of destabilization at cloud base {—d\aN'^/dt) in units of hr~^. Cases 

547 with a Spreading Number S that is much greater than one are indicated in 

548 bold. \M 

549 4 Evaporation rate, in units of h~^, defined here as the negative of the logarith- 

550 mic rate of mass change in the lower depth h of the cloud nto the cloud within 

561 the initial 180 s of simulation in units of s for cases where E > 1 and S* < 1 EHl 

552 5 Dominant modes of evolution observed in simulations. Cases where the Spread- 

553 ing Number S and Evaporation Number E are much greater than one are 

564 indicated in bold and italics, respectively. |3T] 
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Table 1. Spreading number S = — = 





L=100m 


1km 


10km 


?,=0.01g/kg 


1.1x10-^ 


1.1x10-=^ 


0.011 


O.lg/kg 


3.3x10-^ 


0.033 


0.33 


Ig/kg 


13 


130 


1300 
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Table 2. Evaporation number E = = ^ 





L=100m 


1km 


10km 


g,=0.01g/kg 


150 


150 


150 


O.lg/kg 


3.7 


3.7 


3.7 


Ig/kg 


0.037 


0.037 


0.037 
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Table 3. Rate of destabilization at cloud base {—dlnN'^/dt) in units of hr ^. Cases with 
a Spreading Number S that is much greater than one arc indicated in bold. 





L=100m 


1km 


10km 


g,=0.01g/kg 


0.12 


0.17 


0.15 


O.lg/kg 


1.32 


0.83 


1.94 


Ig/kg 


4.00 


11.42 


23.88 
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Table 4. Evaporation rate, in units of h^^, defined here as the negative of the logarithmic 
rate of mass change in the lower depth h of the cloud nto the cloud within the initial 180 s 
of simulation in units of s for cases where E > 1 and < 1 





E = 150 


E = 3.7 


L = 100 m 


5.8 


0.79 


1 km 


4.0 


0.72 


10 km 


1.5 


0.68 



30 



Table 5. Dominant modes of evolution observed in simulations. Cases where the Spreading 
Number S and Evaporation Number E are much greater than one are indicated in bold and 
italics, respectively. 





L=100m 


1km 


10km 


g,=0.01g/kg 


evaporation 


evaporation 


evaporation 


O.lg/kg 


lofting 


lofting 


mixing 


Ig/kg 


mixing 


mixing 


mixing 
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555 List of Figures 

566 1 Radiative energy is transferred from the lower troposphere to the base of a 

567 gravitationally stratified cloud with initial width L, due to a radiative tem- 

568 perature difference AT and deposited into a layer of characteristic depth h at 

569 the base of the cloud. The layer is initially at equilibrium buoyant potential 

560 with respect surrounding clear air at the same level, and it is perturbed 

561 from equilibrium through the deposition of radiative energy into a well-mixed 

562 layer of depth 6z. EH 

563 2 A Schematic diagram of the thermodynamic evolution of a cirrus cloud in 

564 response to radiative diabatic heating. Potential energy flows from the warmer 
566 lower troposphere into the cooler cloud base (red arrow, left). The potential 

566 difference between the cloud and the ground is Ae ~ ^aT^AT where AT 

567 is the effective brightness temperature difference between the cloud and the 

568 lower troposphere. This flow of radiative potential energy perturbs the cloud 

569 from gravitational equilibrium at a rate ^ (red arrow, right). The cloud 

570 acts to restore gravitational equilibrium at rate aAe, resulting in horizontal 

571 spreading of the cloud (blue arrows) . |35] 

572 3 3D surface of cloud as initialized for all simulations ESI 

573 4 Heating rate profiles for = 0.01 g kg~^, = 0.1 g kg~\ and = 1 g kg~^. 

574 Cloud vertical boundaries are marked with a dashed line. 123 
576 5 Cross section of 6^ contours through a cloud with L=lkm and gi=0.1 g kg~^ 

576 [S = 0.033) after s (thin) and 3600 s (thick) of simulation. The initial cloud 

577 boundaries are indicated by the shaded region. |38] 
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578 6 as in Figure [5| but a 3D surface plot of qi after 3600 s of simulation (top), and 

579 a plot of a qi cross section of the cloud (bottom). The initial position of the 

580 cloud is shown in black, while the state of the cloud after 3600 s is shown in 

581 color, with the value of qi denoted by the color scale. Note the rise of cloud 

582 base and the horizontal spreading. 133 

583 7 6^ profile of a cloud with L=10 km and qi=l g kg~^ after 3600 s of simulation. 

584 The initial profile is plotted in a dashed line with horizontal dashed lines 

585 indicating initial cloud base and cloud top. The 6^ profile is calculated as a 

586 horizontal average of all 6^ profiles within an annular region of the cloud. The 

587 inner edge of the annulus is 7.5 km from cloud center and the outer edge of 

588 the annulus is at 9 km from cloud center. HQ] 

589 8 As in Figure [7] but a cross section of 6e contours through a cloud after 

590 s (thin) and 3600 s (thick) of simulation. The initial cloud boundaries are 

591 indicated by the shaded region. Hi] 

592 9 As in Figure [7| but a 3D surface plot of qi after 3600 s of simulation (top), 

593 and a plot of a qi cross section of the cloud (bottom). The initial position of 

594 the cloud is shown in black, while the state of the cloud after 3600 s is shown 

595 in color, with the value of q^ denoted by the color scale. Note that cloud 

596 base remains at roughly the same level and that the cloud bends upward as 

597 it spreads outward. 112] 

598 1 The time evolution of the Spreading Number S and the mixed-layer depth 

599 6z of a cloud with ice water mixint ratio of 1 g kg~^ and width L of 10 

600 km from to 3600 s. The dashed lines indicates slopes of 1/2 and -3 on the 

601 log-log plot as indicated. US] 
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Fig. 1. Radiative energy is transferred from the lower troposphere to the base of a gravita^ 
tionally stratified cloud with initial width L, due to a radiative temperature difference AT 
and deposited into a layer of characteristic depth h at the base of the cloud. The layer is 
initially at equilibrium buoyant potential eeq with respect surrounding clear air at the same 
level, and it is perturbed from equilibrium through the deposition of radiative energy into a 
well-mixed layer of depth 5z. 
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Fig. 2. A Schematic diagram of the thermodynamic evolution of a cirrus cloud in response 
to radiative diabatic heating. Potential energy flows from the warmer lower troposphere into 
the cooler cloud base (red arrow, left). The potential difference between the cloud and the 
ground is Ae ^ ^aT^AT where AT is the effective brightness temperature difference between 
the cloud and the lower troposphere. This flow of radiative potential energy perturbs the 
cloud from gravitational equilibrium at a rate (red arrow, right). The cloud acts to 
restore gravitational equilibrium at rate aAe, resulting in horizontal spreading of the cloud 
(blue arrows). 
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Fig. 3. 3D surface of cloud as initialized for all simulations 
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Heating rate profile, q. = 0.01 g/kg 
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Fig. 4. Heating rate profiles for — 0.01 g kg ^, = 0.1 g kg \ and = 1 g kg ^. Cloud 
vertical boundaries are marked with a dashed line. 
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Fig. 5. Cross section of 9e contours through a cloud with L=lkm and qi=0.1 g kg~^ 
{S — 0.033) after s (thin) and 3600 s (thick) of simulation. The initial cloud boundaries 
are indicated by the shaded region. 
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Fig. 6. as in Figure [5| but a 3D surface plot of qt after 3600 s of simulation (top), and a 
plot of a Qi cross section of the cloud (bottom) . The initial position of the cloud is shown in 
black, while the state of the cloud after 3600 s is shown in color, with the value of qi denoted 
by the color scale. Note the rise of cloud base and the horizontal spreading. 
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Fig. 7. 9v profile of a cloud with L—10 km and Qi—l g kg~^ after 3600 s of simulation. The 
initial profile is plotted in a dashed line with horizontal dashed lines indicating initial cloud 
base and cloud top. The 6^ profile is calculated as a horizontal average of all 6^ profiles 
within an annular region of the cloud. The inner edge of the annulus is 7.5 km from cloud 
center and the outer edge of the annulus is at 9 km from cloud center. 
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Fig. 8. As in Figure a cross section of 6e contours through a cloud after s (thin) 

and 3600 s (thick) of simulation. The initial cloud boundaries are indicated by the shaded 



region. 
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Fig. 9. As in Figure [7| but a 3D surface plot of qi after 3600 s of simulation (top), and a 
plot of a Qi cross section of the cloud (bottom) . The initial position of the cloud is shown in 
black, while the state of the cloud after 3600 s is shown in color, with the value of denoted 
by the color scale. Note that cloud base remains at roughly the same level and that the 
cloud bends upward as it spreads outward. 
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Fig. 10. The time evolution of the Spreading Number S and the mixed-layer depth 6z of a 
cloud with ice water mixint ratio of 1 g kg~^ and width L of 10 km from to 3600 s. The 
dashed hnes indicates slopes of 1/2 and -3 on the log- log plot as indicated. 
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